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Abstract 


We present standard model predictions for the fonr-fermion reaction being one 

of the best detection channels of a low mass Higgs boson prodnced throngh the Higgsstrahlnng 
mechanism at a linear collider. We inclnde leading virtnal and real QED corrections dne to 
initial state radiation and a modification of the Higgs-66 Ynkawa conpling, cansed by the 
rnnning of the 6-qnark mass, for > iJ,'^^~bb. The complete 0{a) electroweak corrections 

to the Z-Higgs prodnction and to the Z boson decay width, as well as remaining QCD and 
EW corrections to the Higgs decay width, as can be calculated with a program HDECAY, are 
taken into account in the double pole approximation. 
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1 INTRODUCTION 


Although the Higgs boson has not yet been discovered, its mass mn can be constrained in the 
framework of the standard model (SM) by the virtual effects it has on precision electroweak 
(EW) observables. Recent global fits to all precision EW data [T] give a central value of 
niH = II 4 I 40 GeV and an upper limit of 241 GeV, both at 95% CL, in agreement with 
combined results on the direct searches for the Higgs boson at LEP that lead to a lower limit 
of 114.4 GeV at 95% CL j2]. These constraints indicate the mass range where the SM Higgs 
boson should be searched for. If the Higgs boson exists, it is most probably to be discovered at 
the Large Hadron Collider, but its properties can be best investigated in a clean experimental 
environment of collisions at a future International Linear Collider (ILC) jS]. 

Main mechanisms of the SM Higgs boson production at the ILC are the Higgsstrahlung reac¬ 
tion 


e+e" ^ ZH, (1) 

the WW fusion 

e+e- ^ UePeW^W* VePeH, (2) 

and the ZZ fusion process 

e+e" ^ e+e-Z*Z* e+e-R. (3) 

The Feynman diagrams of reactions o, © and 0 are depicted in Fig. ^ and lit, 
respectively. 



Figure 1: Feynman diagrams of reactions o, m and ©, respectively. 

The cross section of reaction m decreases according to the 1/s scaling law, while that of 
reaction © grows as ln(s/m|^). Hence, while the Higgsstrahlung dominates the Higgs boson 
production at low energies, the WW fusion process overtakes it at higher energies. The Higgs 
boson production rate through the ZZ fusion mechanism © is by an order of magnitude 
smaller than that of process ©. The production of the SM Higgs boson in the intermediate 
mass range at high-energy colliders through reaction o was already thoroughly studied 
in the literature j^. In the present paper, we contribute further to the theoretical analysis of 
the Higgsstrahlung reaction by taking into account decays of the Z and Higgs bosons, includ¬ 
ing the most relevant EW radiative effects both to the production and decay subprocesses. 
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Preliminary results of the present study have been already presented in our conference paper 

0 


If the Higgs boson has mass in the lower part of the range indicated above, say mu < 140 GeV, 
it would decay dominantly into a bb quark pair. As the Z boson of reaction (HD decays into 
a fermion-antifermion pair too, one actually observes the Higgsstrahlung through reactions 
with four fermions in the hnal state. To be more specific, in the following we will concentrate 
on a four-fermion channel 

e+(pi) + e~{p 2 ) + /i“(P4) + b{p5) + b{pG), (4) 

one of the most relevant for detection of (HD at the ILC, and the corresponding bremsstrahlung 
reaction 


e+(pi) + e (P 2 ) ^ +/^ (P 4 )+ &(P5)+ &(P6)+ 7 (^ 7 )- (5) 

Both in (HD and (ED, the particle four momenta have been indicated in parentheses. While 
hnal states of reactions (HD and © are practically not detectable in hadronic collisions because 
of the overwhelming QCD background, at the ILC, they should leave a particularly clear 
signature in a detector. To the lowest order of the SM, in the unitary gauge and neglecting 
the Higgs boson coupling to the electron, reactions 0 and © receive contributions from 34 
and 236 Feynman diagrams, respectively. Typical examples of Feynman diagrams of reaction 
0 are depicted in Fig. El The Higgsstrahlung ‘signal’ diagram is shown in Fig. la, while the 
diagrams in Figs, lb and Ic represent typical ‘background’ diagrams. The Feynman diagrams 
of reaction © are obtained from those of reaction © by attaching an external photon line 
to each electrically charged particle line. Reactions 0 and © are typical examples of the 
neutral current four-fermion reactions and the corresponding hard bremsstrahlung reactions 
which had been already studied for LEP2 and TESLA TDR j3] in the case of massless fermions 
P. P and without neglecting fermion masses P. p. Reaction © was in particular studied 
in ng. It was also considered, among other four-fermion reactions relevant for the Higgs 
boson production and decay, in a study of the Higgs boson searches at LEP2 performed in 
HU, where references to other works dedicated to the subject can be found, too. 



Figure 2: Examples of Feynman diagrams of reaction 0: (a) the double resonance ‘signal’, 
(b) and (c) ‘background’ diagrams. 

In order to match the precision of data from a high luminosity linear collider, expected to 
be better than 1%, it is necessary to include radiative corrections in the SM predictions for 
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reaction ©• Calculation of the complete electroweak (EW) 0{a) corrections to a four-fermion 
reaction like 0 is a challenging task. Problems encountered in the hrst attempt of such a 
complete calculation were described in [T2j. Substantial progress in a full one-loop calculation 
for —*• udfi~u^ was reported by the GRACE/l-LOOP team [T2] and recently results of the 

hrst calculation of the complete EW 0(a) corrections to charged-current > 4 fermion 

processes has been presented ^1]. Another step towards obtaining high precision predictions 
for the Higgs boson production, accomplished in ra. was calculations of EW corrections 
to —>• uuH, a process related to 0. At the moment, however, there is no complete 

calculation of the corrections to neutral-current e''“e“ —> 4 fermion processes available that 
could be included into a Monte Carlo (MC) generator. Therefore it seems natural to include 
the leading QED corrections and to apply the double pole approximation (DPA) for the EW 
0(a) corrections, as it had been done before in in the case of W-pair production at LEP2. 
This means that we will include the so called factorizable EW corrections to & ZH production 
o and to subprocesses of the Z and Higgs decay, which are available in the literature. 

To the lowest order of the SM, the cross sections of reactions 0 and © can be computed with 
a program ee4fg na. On the basis of ee4fg, we have written a dedicated program eezh4f 
that includes the following radiative corrections to 


• the virtual and real soft photon QED corrections to the on-shell ZH production O, a 
universal part of which is utilized for all Feynman diagrams of reaction (ji} and combined 
with the initial state hard bremsstrahlung correction of 0, which we refer to as the 
initial state radiation (ISR) QED correction; 

• the weak corrections to the Z-Higgs production |THI, he. the infrared (IR) hnite 

part of the complete 0{a) EW correction to 0 and full EW corrections, including hard 
photon emission, to the Z decay width 1213], [23 DPA; 

• the correction Higgs boson decay width, as can be computed with HDECAY 


The correction to the partial Higgs boson decay width T n^bb is split into two parts: one part 
is the tree level partial width parametrized in terms of the running 6-quark mass and 

the other part is the correction term that contains the remaining QCD corrections and the 
bulk of EW corrections, as described in |22|- The hrst part is included in the calculation of 
reactions 0 and through the corresponding modihcation of the H —>■ bb Yukawa coupling, 
while the second part is taken into account in the DPA. 

The non-factorizable corrections to 0 have not been included in eezh4f yet. Although, a 
comparison with the related process of ZZ production and decay, where the non-factorizable 
corrections largely cancel each other in the total cross section |22| , [211 , may suggest that they 
are also numerically suppressed for the ZiY-mediated four-fermion hnal states if decay angles 
are integrated over, they should be included in the program if one wants to study differential 
cross sections, or if one wants to improve predictions for invariant-mass distributions and 
similar quantities. 
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2 SCHEME OF THE CALCULATION 


The corrections listed in the introdnction are taken into acconnt in the radiatively corrected 
total cross section of 0 according to the master formnla 


dcr / d(JBorn+r.m. T 


/ Ey<Ecut 


1 QED ISR I 
virt+soft,univ. 


' Ey>Ecut 


1 QED ISR 
•^^hard 


I f J^EW 

+ J 


EW DPA 
finite * 


By d(TBorn+r.m. we denote the effective Born approximation 


dcr 


Born+r.m. 


1 


M 


■4/ 


Born 


+ 2Re{MZlSM^L)]d^4f, 


( 6 ) 

(7) 


where is the matrix element of reaction (j3]) obtained with the complete set of the lowest 

order Feynman diagrams, is the correction to the amplitnde of the ‘signal’ diagram of 

Fig. Ha) dne to modihcation of the lowest order Higgs -66 Yukawa coupling caused 

by the running of the 6 -quark mass, d$ 4 / is the Lorentz invariant four-particle phase space 
element and 4/ is a shorthand notation for the hnal state of reaction 0. The correction 
can be written as 


^fY0) , AO) \ r ^-9^ +PuP3aMI 

= V Pi [A 7/. + 71^ '7m75 ) U ps ^ 

^ ^ Dz{P34)Dh{P56) 

X u (pa) {ai-f,, + a2-ful5) V (ps) 6^^“^ h(p6)n(p5). 


( 8 ) 


where and are the lowest order invariant amplitudes of the on-shell Higgsstrahlung 
reaction 0 HHI which read 


( 0 ) 

4 ( 0 ) 9hzz 

af\ z = l,2 

(9) 

s-Ml 

and Dz{p 3 a) and Dh^p^q) are denominators 

of the Z and Higgs boson propagators 


Dz{p3i) = (P3 + P4f - M|, 

Dh{P5g) = {Pb+Pef - Mh, 

(10) 

with the complex mass parameters 



Ml = ml - imzTz, 

Mjj = 

(11) 


The ‘hxed’ total widths Tz and Ti;/ of the Z and Higgs bosons have been introduced in order 
to avoid singularities in resonant regions of the corresponding propagators. We have used 
the complex mass substitution also in Eq. Q, although the propagator factor of Eq. 0 can 
never become resonant as we want to keep open a possibility of performing calculations in the 
complex mass scheme j7j, which preserves Ward identities. 


To lowest order of the SM, the vector and axial-vector couplings of the Z boson to leptons, 
af\ a 2 °\ the Higgs boson coupling to Z, g^^z of Fq. Q and the Higgs -66 Yukawa coupling, 
gmb, are given by 


a® 


4g^ - 1 

4:S\yC]V 


ew: 


= 


evE 


dSiyCiy ’ 


(0) 

9hzz 


ewmz 

SwCw 


(0) 

Qmb 


(^wmb 
2swmw ’ 
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with the effective electric charge ew = (47ravv)^'^^ and electroweak mixing angle 6w defined by 

aw = \/2G^rri^s‘^/TT, and = 1 - (13) 

where niw and mz are the physical masses of the W and Z boson, respectively. This choice, 
which is exactly eqnivalent to the G^-scheme of HHI, is in the program referred to as the 
‘fixed width scheme’ (FWS). We have introduced usual shorthand notation sw = sin6*w and 
cw = cos 9w in Eqs. m and m- 

The total widths of Eq. oa are calculated numerically in the framework of the SM: T^ is 
computed with a program based on reference HP and Th is obtained with a program HDECAY 
ra, both programs including radiative corrections. In order not to violate unitarity it is 
important that we include exactly the same corrections in the amplitudes of the partial decay 
widths Tz^f,+i,- and 

The program HDECAY includes the full massive NLO QCD corrections close to the thresholds 
and the massless corrections far above the thresholds. Both regions are related with a 

simple linear interpolation equation which, for the partial Higgs decay width T reads 

= (l - ^(mb ), with r = (14) 

where large logarithms are resummed in the running 6-quark quark mass in the MS renormal¬ 
ization scheme Wib{mH)- 


By inserting the representations 

T H^bbimbimu)) = V^^^^iimbimH)) + H^bb 

T H^bb {mb) = T^H^bb {^b) + at H^bb {mb) , 
with the lowest order partial width into 66-quark-pair given by 


T 


( 0 ) 

H^bb 


{mb) = 3 


( 0 ) 2 
Snbb 

Svr 


mH\l- 


4m 


2\ i 


ml 


into Eq. (HI, we obtain the representation 


^ H^bb — ^'^mlbb + ^^H^bb 


for T H^bb- The hrst term on the right hand side of Eq. (j18jl 


mn 


(15) 

(16) 


(17) 


(18) 


rSiiS = (l - + r^rgtiS (rnt) (19) 

contains the correction due to running of the 6-quark mass, while the second term, 
incorporates the remaining QCD corrections and EW corrections that are taken into account 
in HDECAY. 


Now, if we write the radiatively corrected Higgs-66 Yukawa coupling in the form 

hBbb = 9 nib T ^d^Hbi + ^9Hbb 


( 20 ) 
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and the corrected partial Higgs decay width as 


p _ _ p(o) 

H^bb — 


H^bb 




1 + 2 Re 


^9 mi + ^dHbb 

^Hbb 


( 21 ) 


then we will obtain the following expressions for the radiative corrections and 5gmh 

assuming that they real, 


89 


r.m. _ 

Hbb ~ 


„(«> 

2rSUm) 


^9Hbb 


■p _ pr.m. 

1 H^bb - >- H^bb 

2rgU(mi) 


^Hbb' 


( 22 ) 


We see that taking into account the correction in Eq. dZD is actually equivalent to 

replacing mb in the lowest order Higgs-66 Yukawa coupling of Eq. m by some effective value 
of the 6-quark mass. The same modihcation of the Higgs-66 Yukawa coupling is done in the 
soft and hard bremsstrahlung corrections represented by the second and third term on the 
right hand side of Eq. ©, respectively. 


After having arranged the radiative corrections in the manner just described, the soft brems¬ 
strahlung contribution of Eq. can be written as 


da5f,° 


ISR 


virt+soft,univ. 


= da 


^QED 


ISR 


Born+r.m.'-'virt+soft,univ. 5 


(23) 


where the correction factor 


^QED ISR 
'^virt+soft,univ. 


27r2 


, mi 


, 2Ecut 3 s 
In —^ -h - In —^ 




mi 


(24) 


combines the universal IR singular part of the 0{a) virtual QED correction to the on-shell 
ZH production process (CD with the soft bremsstrahlung correction to 0, integrated up to 
the soft photon energy cut Ecnt- In Eq. fj2d|l . e is the electric charge that is given in terms of 
the fine structure constant in the Thomson limit ckq, e = (dvrcxo)^^^. 


The third term on the right hand side of Eq. (jHD represents the initial state real hard photon 
correction to reaction 0, i.e. the lowest order cross section of reaction 0 with the photon 
energy cut > E^nt, which is calculated taking into account the photon emission from the 
initial state particles. It has been checked numerically that the dependence on Ecnt cancels in 
the sum of the second and third term on the right hand side of Eq. 0, provided that the real 
photon coupling to the initial state fermions is parametrized in terms of ckq, too. 

Finally, the last integrand on the right hand side of Eq. (jHD, is the IR hnite part of 

the virtual EW 0(a) correction to reaction in the DPA. It can be written in the following 
way 


1 EW, DPA 
finite 


2s 




r 


non—univ. 


QED 


+ 2Re ( 


dpa^Mdpa) I 


d$ 


4 /) 


(25) 


where the lowest order matrix element and the one-loop correction SMppA in the DPA, 

which are given below, are calculated with the projected four momenta ki, i = 3,..., 6, of the 
hnal state particles, except for denominators of the Z and Higgs boson propagators. The 
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projected four momenta are obtained from the four momenta pi, i = 3, of reaction 0 
with a, to some extent arbitrary, projection procedure which will be dehned later; 
denotes the IR finite non universal constant part of the 0{a) QED correction that has not 
been taken into account in Eq. (EH). It reads 

^QED ISR / TT^ \ 

Rnon-univ. = ^ I Y “ H ' ^^6) 


The matrix elements and SMppA for the Z and Higgs boson production and decay of 

Eq. are given by 


Mdpa = -V (Pi) (dl®7/. + u (pa) 


6 Mdpa — 


V (pi 

+ (H 41 PI;, + H42P2^) jtz'y5 

+ V (pi) + ^^7^75) u (pa) (sN^ Nh + SN, 


Nh, 

'z 

(27) 

U {P2) Np 

(28) 


Dz{pu)DH{pm) ’ 


with denominators of the Z and Higgs boson propagators defined by Eq. (uni) and amplitudes 
of decay subprocesses given by 

N^ = u (ki) (af V (k^) , Np = gplb “ (^e) v ih) (29) 

in the lowest order and 

= u (/C 4 ) (hai 7 ^ + ^ 027 ^ 5 ) V {ks), 5Np = Sgpbb u {k^) v {kb) (30) 


in the one-loop and higher order. Let us note that the invariant amplitudes hui, 6 a 2 and 
Sgpbb contain the QED correction factors, which include the hard photon emission. We have 
neglected the so-called longitudinal part of the Z-boson propagator proportional to P 34 P 34 
in Eqs. (EH) and (EED- Keeping this part would lead to some ambiguity, as the propagator 
couples to the hnal state currents and 6 N^ which are parametized in terms of the projected 
momenta k^ and k^. The neglected term is of the order of m^, which is consistent with neglect 
of the 0{m^) terms in calculation of the one-loop EW decay amplitudes 6 ai and 6 a 2 of Eq. (IHIH) . 

The invariant amplitudes of Eq. EH|)^ 

= 5Hj(s, cos^), Aji = Aji{s, cos 9), i = l,2, j = 3,4, (31) 

represent the IR finite parts of the EW one-loop correction to the on-shell Z-Higgs production 
process O- They are complex functions of s = (pi + ^ 2 )^ and cos^, 9 being the Higgs boson 
production angle with respect to the initial positron beam in the centre of mass system (CMS), 
computed with the program eezh4f. The latter makes use of the program worked out in EHl 
and the package FF 2.0. The package FF written by G. J. van Oldenborgh allows to evaluate 
one-loop integrals |2S|. Note that we have changed a little the notation in Eq. (EHl) with respect 


^Only amplitudes which do not vanish in the limit rrie —> 0 are included. 








to that of Eq. (2.5) of ^H]- The IR hnite contributions to the EW one-loop form factors of the 
Zll-veitex Sai, 6ai and the correction Sgnbb of the Hbb-yertex of Eq. (EDD are also calculated 
numerically following references |2I] and [221 • 

As the computation of the one-loop electroweak amplitudes of Eq. m slows down the MC 
integration substantially, a simple interpolation routine has been written that samples the 
amplitudes at a few hundred values of cos 6 and then the amplitudes for all intermediate 
values of cos6' are obtained by a linear interpolation. This gives a tremendous gain in the 
speed of computation, while there is practically no difference between the results obtained 
with the interpolation routine and without it. 


We end this section with a description of the projection procedure we have applied in the 
program. The projected momenta kz and fc,, of Eqs. (121, (El and El are obtained from 
the four momenta pi, i = 3,..., 6, of the hnal state fermions of reaction o with the following 
projection procedure. 


First the on-shell momenta and energies of the Higgs and Z boson in the CMS are hxed by 


kn — 
Eh = 



kn 

kz 


r P5 + Pe 
kn -p—— 

\P5 +Pel 

—kn, Ez = a/s — Eh- 


(32) 


Next the four momenta p and p^ of reaction © are boosted to the rest frame of the bb- and 
/i+/i“-pairs, respectively, where they are denoted by p'^ and p'^. The projected four momenta 
k'-, i = 5, 6, in the rest frame of the 66-pair are then obtained by the kinematical relations 


k'r> 

k\ 


A2 


2 mH 


— —k'^i 




E’i = 


k'r 


P'b 


P'b 


1 

,2 2 

k'i +mi] , i = 5, 6. 


(33) 


Similarly, one obtains /cg and k'^ in the rest frame of the pZp -pair using 


A 2 (m|, m^) 


k’. = 


p's 


2mz 


p's 


kU = -k's, 


E'i = {k\+mi] , i = 3,4. 


(34) 


The four momenta k[, i = 3 ,..., 6, are then boosted back to the CMS giving the projected four 
momenta ki, i = 3,..., 6 of p~^, p~,b and 6 that satisfy the necessary on-shell relations 

kl = kl = ml, kl = kl = ml, 

{ks + kif = ml, (fcs -F kef = m]j. 

The actual value of cos 6* in Eq. EH) is given by cos^ = k^^/\kH\- The described projection 
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procedure is not unique. As the Higgs boson width is small, the ambiguity between different 
possible projections is mainly related to the off-shellness of the Z boson and is of the order of 
aVzl{T^mz)- 


3 NUMERICAL RESULTS 


In this section, we will present a sample of numerical results for reaction m which have been 
obtained with the current version of a program eezh4f. Computations have been performed in 
the FWS with the Z boson mass, Fermi coupling and fine structure constant in the Thomson 
limit as the initial SM EW physical parameters j2tij : 

mz = 91.1876 GeV, = 1.16639 x lO"*^ GeV^ oq = 1/137.03599976. (35) 

The external particle masses of reaction 0 are the following: 

me = 0.510998902 MeV, = 105.658357 MeV, mb = 4.4 GeV. (36) 

For definiteness, we give also values of other fermion masses used in the computation: 

mr = 1.77699 GeV, = 1.5 GeV, mt = 177.7 GeV, 

mu = 62 MeV, md = 83 MeV, = 215 MeV. (37) 


The light quark masses of Eq. (EZI, together with = 0.123, reproduce the hadronic contri¬ 
bution to the running of the fine structure constant. 

Assuming a specific value of the Higgs boson mass, the W boson mass and the total Z boson 
width are calculated in a subroutine based on izu. while the total Higgs boson width is cal¬ 
culated with HDECAY Q. We obtain the following values for them for m^ = 115 GeV and 
parameters specihed in Eqs. (jSSHSIl) 


mw = 80.40844 GeV, Tz = 2.50393 GeV F^ = 2.8542 MeV. (38) 


The total ‘signal’ cross section of 0 in the narrow width approximation is plotted on the left 
hand side of Fig. IHlas a function of the GMS energy. The solid curve shows the Born cross 
section 

■p(O) -pio) 

NWA _ (0) ^ ^ H^bb 

^Born ^ e+e~^ZH p p 

t Z t H 

The dashed curve shows the cross section including the correction due to running of the 6-quark 
mass 



pr.m. 

NWA ^ (0) ^ ^ H^bb 

^ Born+running mj, ^ e+e~^ZH p p 

t Z t H 


(40) 
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Vs (GeV) 



Vi (GeV) 


Figure 3: The ‘signal’ total cross section of reaction © in the NWA as a function of the CMS 
energy. The plots on the left show the Born cross section of Eq. (IHHl) (solid line), the cross 
section of Eq. (Unj) (dashed line) and the cross section including all corrections as given by 
Eq. (11^ (dotted line). The plot on the right shows the full relative correction with respect to 
the Born cross section. 


Finally, the dotted curve shows the cross section including the complete EW corrections to the 
Z-Higgs production and to the Z boson decay width, as well as the QCD and EW corrections 
to the Higgs boson decay width obtained with HDECAY and the ISR. The latter has been taken 
into account in the structure function approach by folding the cross section 


dcr 


NWA 


= dcr, 


e+e-^ZH- 




H^bb 


(41) 


Tz Th 

with the structure function Fg/' [x, Q^) and integrating it over the full angular range according 
to 

d^NWA 


a 


NWA 

Born+all correc. 


= / dxi 


pi i\WA 

dxaF^g^ F“ (x2,(5^) j^^dcose^-^ {xipi,X2P2), (42) 


where xipi {X 2 P 2 ) is the four momentum of the positron (electron) after emission of a collinear 
photon. The structure function F£f (x, Q^) is given by Eq. (67) of j2ZI, with ‘BETA’ chosen 
for non-leading terms. The splitting scale Q^, which is not hxed in the LL approximation is 
chosen to be equal s = (pi +^ 2 )^- The corresponding complete relative correction is plotted 
on the right hand side of Fig. El It is large and dominated by the correction to the Higgs-66 
Yukawa coupling caused by the running of rrih. 


How the ‘background’ Feynman diagram contribution and off-shell effects change the Hig- 
gsstrahlung signal cross section is illustrated in Fig. 01 where on the left (right) hand side we 
plot the total Born cross section of reaction © as a function of the CMS energy in the energy 
range from 0.12 - 1 TeV (0.19 - 0.3 TeV) for a few different values of mn- The Z-Higgs pro¬ 
duction signal that is clearly visible in Fig. 01 as a rise in the plots above the solid line, which 
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represents the cross section of reaction dl without the Higgs boson exchange contribution, 
decreases while the Higgs mass is growing. For mn = 150 GeV, it becomes already rather 
small, while it is hardly visible for tuh = 160 GeV. This effect is caused by the substantial 
grow of the total Higgs boson width F/^ with the grow of the Higgs boson mass rriH- The 
latter results in a decrease of the branching ratio Vu, as with growing mn the Higgs 
decay into hF-boson pair starts to dominate. The first bump in the plots on the left hand side 
of Fig. m reflects the double Z production resonance. 



Figure 4: The total Born cross section of reaction 0 as a function of the GMS energy in the 
energy range from 120 - 1000 GeV (left) and from 190 - 300 GeV (right) for different values 
of TUh- 


The total cross sections of reaction Q including radiative corrections are plotted in Fig. El 
as functions of the GMS energy in the range from 190 - 300 GeV for a few selected values 
of the Higgs boson mass. The plots on the left include the correction to the Higgs-6& 

Yukawa coupling caused by the running of the Yquark mass, he., they depict the first integral 
on the right hand side of Eq. 0. The correction is large and it reduces the Higgsstrahlung 
signal substantially as compared to the plots on the right hand sight of Fig. E) The plots on 
the right hand side of Fig. El include all the corrections as dehned in Eq. (jH)). As compared to 
the corresponding plots on the left, they are shifted upwards with respect to the solid curve 
representing the Born cross section of reaction 0 without the Higgs boson contribution. This 
is caused by taking into account the contribution from the initial state bremsstrahlung with 
the hard photon momentum integrated over the full phase space. As expected, inclusion of 
the ISR smears the Z-Higgs production signal. 

The size of EW corrections corresponding to the fourth term on the right hand side of Eq. © 
can be read from Fig. El where the total cross section of reaction Q for mu = 115 GeV 
(left) and niH = 130 GeV (right) including all the corrections (dashed-dotted line) and the 
cross section including the QED ISR corrections and the correction due to the running of 
mb, corresponding to the first 3 terms on the right hand side of Eq. (jH]), (short dashed line) 
are compared. The plots of the Born cross section of reaction 0 without the Higgs boson 
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Figure 5: The total cross section of reaction Q as a function of the CMS energy for different 
values of mn- The plots on the left include the correction to the Higgs-66 Yukawa 

coupling caused by the running of the 6-quark mass. The plots on the right include all the 
corrections of Eq. 



Figure 6: The total cross section of reaction (jU for niH = 115 GeV (left) and mu = 130 GeV 
(right) including all the corrections (dashed-dotted line) as compared to the cross section 
including the QED ISR corrections and the correction due to the running of mb, corresponding 
to the hrst 3 terms on the right hand side of Eq. (jH)), (short dashed line). The plots of the 
Born cross section of reaction dH) without the Higgs boson exchange (solid line) and the cross 
section of 0 including solely the correction due to the running of mb (dotted line) are shown, 
too. 
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exchange (solid line) and the cross section of (jl)) inclnding solely the correction due to the 
running of rrih are shown, too. 


4 SUMMARY AND OUTLOOK 


We have presented the SM predictions for a four-fermion reaction which is one of the best 
detection channels of a low mass Higgs boson produced through the Higgsstrahlung mechanism 
at a linear collider. We have included leading virtual and real QED ISR corrections to all the 
lowest order Feynman diagrams of reaction 0. We have modihed the lowest order Higgs- 
66 Yukawa coupling in reactions © and © by including in it the correction to the Higgs 
boson decay width caused by the running of the 6-quark mass, which can be calculated with 
a program HDECAY. The complete 0{a) EW corrections to the Z-Higgs production and to the 
Z decay width, as well as the EW and remaining QCD corrections of HDECAY have been taken 
into account in the double pole approximation. Moreover, we have illustrated how the Higgs 
boson production signal is visible in the CMS dependence of the total cross section of reaction 
© for low Higgs masses and how it almost disappears for the Higgs mass approaching the 
lE-pair decay threshold. 

Some further work is required in order to include the non-factorizable corrections to reaction 
©. Although, by comparison to a related process of the ZZ production [221, might 

expect that the non-factorizable corrections may largely cancel each other in the total cross 
section and may be numerically suppressed in differential cross sections integrated over the 
decay angles also for reaction ©, they should be included in the program in order to properly 
take into account radiative corrections due to hnal state real photon radiation which, in the 
present work, have been included inclusively only in the EW one-loop corrections to the Z 
and Higgs boson decay widths. To avoid possible negative weight events which may occur in 
the MC simulation one should exponentiate the IR sensitive terms. 
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